Nonlinearity from quantum mechanics: 
Dynamically unstable Bose-Einstein condensate in a double-well trap 
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We study theoretically an atomic Bose-Einstein condensate in a double-well trap both quantum 
mechanically and classically under conditions such that in the classical model an unstable equilibrium 
dissolves into large-scale oscillations of the atoms between the potential wells. Quantum mechanics 
alone does not exhibit such nonlinear dynamics, but measurements of the atom numbers in the 
potential wells may nevertheless cause the condensate to behave essentially classically. 
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Even though quantum mechanics is supposedly the 
fundamental theory, a classical description succeeds for 
a large host of systems. When the classical model is 
nonlinear, apparent discrepancies with linear quantum 
mechanics readily arise. We have earlier [l| studied sta- 
tionary states of a Bose-Einstein condensate (BEC) in 
an optical lattice. Classical theory predicts a nonlinear 
soliton, while quantum mechanics does not. Of course, 
solitons are seen experimentally [4]. Our resolution of 
this dilemma is that the measurements probing the soli- 
ton will bring it about [l|. 

Here we study a time-dependent case in which quan- 
tum mechanics and classical mechanics seem to disagree, 
namely, dynamical instability. Classically, the signature 
of an instability is that a small deviation from an unsta- 
ble equilibrium grows exponentially in time. Quantum 
mechanically, under unitary time evolution, the distance 
between any two states remains constant. 

A Bose-Einstein condensate in a symmetric double- 
well potential makes our example. A convenient quan- 
tum description exists [3|, and can be easily solved nu- 
merically for a large number of atoms [J| . There is a cor- 
responding classical model, too, exactly as in quantized 
and classical descriptions of an electromagnetic field. It 
is nonlinear, and has been discussed in particular as an 
example of population trapping [3, [a] , asymmetric oscil- 
lations of the atoms between the two wells. Population 
trapping has also been seen experimentally [6| . The clas- 
sical system may have an unstable equilibrium, and a 
smallest deviation from the equilibrium can lead to large- 
scale oscillations of the atoms between the two sides of 
the potential well [7[. The question is, what does quan- 
tum mechanics have to say about such oscillations. 

We describe a preparation that leaves the classical 
model in an unstable state. The oscillations of the atoms 
differentiate between the two potential wells. Since quan- 
tum mechanics strictly preserves the symmetry between 
the sides of the trap, there seemingly cannot be any such 
oscillations. But the broken symmetry in the classical 
case [l|, Isl-llOt once more [l|, |1, Q correctly suggests that 
including a model for the observations of the atoms will 
cure the discrepancy. We do this much as before Q, 
and solve the theory using quantum trajectory simula- 
tions [ll| - [l4| . The quantum system, thus amended, will 



not only^display the classical behavior, but, in analogy 
to Ref. pj, the measurements are seen to literally cause 
the classical dynamics. 

We study a double-well trap for bosonic atoms 
within the conventional two- mode approximation 13| . The 
Hamiltonian is 

^ = -J{blbl + b\br) + U{b].b].brbr + b\b\blbi) . (1) 

Here 6; and 6^ are the annihilation operators for the 
atoms in the "left" and "right" wells, J > is the am- 
plitude for tunneling between the wells, and U charac- 
terizes the strength of the atom-atom interactions. The 
total number the of atoms N = 6^,6^ + b\bi is a constant 
of the motion, and its value is denoted by A^. 

The corresponding classical field theory is obtained by 
taking the Heisenberg equations of motion for the boson 
operators, e.g., ib^ = — Jbi-\-2Ub\b\br , and declaring that 
in the equations of motion the operators are c- numbers, 
no longer operators. We then have, for instance, ibr = 
— Jbi -\- 2U\br\'^br- This classical system has only two 
relevant dynamical variables [5| that we pick as 

z = {\br\^^\bi\^)/N, ^ = arg(6,6|), (2) 

relative population imbalance and phase difference of the 
condensates between the right and left traps. The vari- 
ables z and Lp make a canonical pair, and their evolution 
is governed by the classical Hamiltonian 
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-2 JVi - 2^ cos ip + x{l + z^)]X = NU . (3) 



The variables z and ^p scale away the dependence on 
total atom number A^, but A^ is still is part of the classical 
Hamiltonian ([3]) in that x ^ ^ ■ If ever, the classical 
approximation can be accurate only for a large number 
of atoms, when the discreteness of the atoms may be 
ignored. In the classical limit of the quantum problem 
embodied in Eq. ^ one should simultaneously take A^ ^ 
1, and scale the atom-atom interaction strength U so that 
X = NU remains constant. As the interaction strength 
U can be adjusted using the Feshbach resonance [l5| and 
the atom number can also be varied experimentally, the 
two-well model represents a realistic system that can be 
tuned between the quantum and classical limits. 
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FIG. 1. (Color online) Phase space portrait for the classical 
model of a two-well system for the parameter values xl J ~ 
— 1.5. The thick line marks the homoclinic orbits; the arrows 
indicate the direction of flow in phase space. 

The classical Hamiltonian ([3|) has two stationary states 
in the region (2, t/j) £ [—1,1] ® [0, 27r), namely, {z,ip) = 
(0, 0) and (0,7r). Both are dynamically stable for \)^l J < 
1; for x/J > 1 the steady state (z, ip) = (0, n) is unstable, 
while for x/J < —1 the unstable state is {z,ip) — (0,0). 
Figure [T] shows a phase-space portrait, contours of con- 
stant T-L^z, if). Since the Hamiltonian is a constant of the 
motion, the phase-space coordinate {z,ip) is constrained 
to move along a constant-energy contour. Figure [T] is 
for the case with x/J = —1.5, whereupon the point 
(z, if) — (0, 0) is a dynamically unstable equilibrium, a 
hyperbolic fixed point. There are two directions in phase 
space in which the system first recedes exponentially in 
time from the point (0, 0) and then returns to the fixed 
point, completing a loop called homoclinic orbit. 

To facilitate an unbiased comparison between classical 
and quantum dynamics, we envisage an experimentally 
feasible preparation that works in the same way for both 
cases. Thus, the system is first prepared to the lowest- 
energy state in the presence of repulsive atom-atom in- 
teractions, x/J > 1, and at the time i = the sign of 
the atom-atom interaction is suddenly fiipped. 

Classically, the lowest-energy state for positive interac- 
tions is (z, if) = (0, 0), with half of the atoms on each side 
of the trap and the same phase for the condensates on 
both sides; for x/J < —1 the same state is the unstable 
equilibrium. Whether in numerical or physical experi- 
ments, a classical system will not stay in a dynamically 
unstable state, but some form of noise will invariably 
launch the instability. In the example of Fig. [T] even a 
minute amount of noise can cause three distinct macro- 
scopic behaviors [3 ■ If the system starts inside the ho- 
moclinic orbit with z > 0, in the absence of further noise 
it will stay inside the same orbit. Correspondingly, the 
atoms stay predominantly in the right trap. The same 
situation but favoring the left trap occurs if the system 
starts inside the homoclinic orbit with z < 0. The third 
alternative is that the system starts (slightly) outside of 
the homoclinic orbits. Then symmetric oscillations of the 
atoms from side to side will result. 



In the quantum case, the preparation puts the system 
in the ground state for repulsive atom- atom interactions. 
This state interpolates from the situation with all atoms 
in the single-particle ground state of the double- well trap, 
IV') oc {bl + bl)^\0), to the number state with half of the 

atoms in each side, \ip) oc hi "bj |0), as the interac- 
tion strength ranges from [/ = to NU/J — !■ 00 0,[J. It 
is not a stationary state after the sign of the interaction 
is flipped. However, since both the ideal initial state and 
the Hamiltonian are invariant under the exchange of the 
site labels I and r, Hamiltonian time evolution alone can- 
not produce any overt difference between the two sides 
of the trap cither: every quantum expectation value re- 
mains unchanged under the exchange of I and r. More- 
over, and in contrast to the classical case, a small error 
in the initial state will not cause runaway oscillations of 
the atoms between the two sides. To find a counterpart 
of classical instability we look elsewhere, into measure- 
ments. 

In our thought experiments we assume in close anal- 
ogy with Refs. Q and [la| that, perhaps by utilizing far- 
off resonant light scattering, a situation has been set up 
whereby two photon detectors, "left" and "right", give 
photon counts as a result of the presence of the atoms 
in the left and right traps at the rates Rr = r{blbr) 
and Ri = r(6j6;). The constant F could depend on ge- 
ometry of the experiment and the tuning and intensity 
of the probe light. Here we make a technical assump- 
tion that removes various secondary complications from 
the picture: no photons are missed, but there is instead 
a one-to-one correspondence between photon detections 
and their back-action on the condensate. Formally, we 
model the photon detections by the operators 
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such that the Liouville-von Neumann equation of motion 
of the density operator of the BEC, p, is amended by an 
addition in the usual Lindblad form [14| . 



P=-JiiH,p] 
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{2UpL\ - l\l,p - pL\L,) . (5) 



The equation of motion ([5]) is unraveled using standard 
quantum trajectory simulations jlll4l4{ . For reasons of 
efficiency and accuracy the practical algorithm is more 
complicated [l^, Hjl than our conceptual discussion [ll[ , 
but the idea is as follows. Suppose we have a state vector 
for the system at time i, \i>{t)). To obtain the state vector 
a short time dt later, we first evolve the state according 
to 
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to obtain \'4>{t + dt)). In our particular example the com- 
bination of the Lindblad operators give ^^ L\Li = ^FiV, 
so that the effect of the term added to the Schrodinger 




equation is simply to damp the state vector by the fac- 
tor e~2 *. At the end of the integration over dt, the 
probabihty for a photon count on the right detector is cal- 
culated, dP,. = dtRr = dtT{b\.hr), and dPi analogously. 
Using a random number generator, the algorithm decides 
whether a photon count happens on either detector. If 
the decision is that the right detector r clicked, the state 
vector at time t + dt will be \il;{t + dt)) = Lr\ip{t)) oc 

if the detector / reported, the result is 
blbi IV'(^)); if neither detector reported. 



\ip{t + dt)) ex 

the state at t + dt is \ip{t + dt)) = \ip{t + dt)); and for a 
short enough time steps dt the possibility that both de- 
tectors clicked is negligible. Finally, the state \il}{t + dt)) 
is normalized, and the next time step commences. 

A collection of stochastic state vectors obtained in this 
way could be used to compute expectation values over the 
density operator p |lll4l4|. More to the present point, it 
can be argued that the sequence of photon counts pro- 
duced along with an individual realization of a quantum 
trajectory \ip{t)) is a representative example of a sequence 
of photon counts that one would see in one run of a con- 
forming experiment. This idea permeates the literature 
on quantum trajectory simulations [lj|, but for complete- 
ness we enunciate the underlying (meta)physical assump- 
tion [l7| . Namely, by construction of the quantum tra- 
jectory algorithm the multitime correlation functions for 
the photon counts from repeated runs of the simulation 
would be the same as the correlation functions from re- 
peated real experiments as predicted from quantum me- 
chanics under the Markov approximation and the quan- 
tum regression theorem |14| . If one reasonably assumes 
that it is not possible to tell the difference between the re- 
alizations of two stochastic processes with the same cor- 
relation functions, quantum trajectory simulations and 
experiments should produce indistinguishable sequences 
of photon counts. 

We may now consider what would happen in an in- 
dividual experiment in which the quantum system starts 
from the state whose classical counterpart is dynamically 
unstable. If we were to analyze an actual experiment, in 
the present setup the only information available would 
be the timing of the photon counts. We would have to 
reconstruct the counting rates at both detectors from the 
observed photon counts, and infer from those the num- 
bers of the atoms in each potential well. In a simula- 
tion we have the shortcut that the instantaneous photon 
counting rates Rr and Ri are on hand, and we may use 
them to find a representation of the population imbal- 
ance. Analogously to Eqs. ([2|) we compute 



zit)^m)\iblbr-blbi)m))/N. 



(7) 



In the simulations we also define the phase difference be- 
tween the condensates by closely mimicking Eq. ([2]), 




^{t) = avg{^{t)\brbl\^{t)) 



(8) 



FIG. 2. (Color online) Time dependence of population im- 
balance z{t) and the phase space trajectory z{ip) from three 
quantum trajectory simulations that differ in the frequency 
of the observations by photon counting; from top to bottom, 
NT/ J = 1, 100, and 10''. The trajectory in phase space is 
plotted as discrete points, 5000 of them; the solid line shows 
the same homoclinic orbits as Fig[l]. 

It is possible to measure experimentally the phase dif- 
ference between the condensates [I^, but in our scheme 
ip(t) from Eq. (jH]) is an auxiliary quantity with no direct 
operational meaning. 

Example results are given in Fig. [51 Here we have 
the atom number N = 10* and the interaction strength 
x/J = NU/J ~ —1.5. Each data set shows the outcome 
of one quantum trajectory simulation. The difference is 
in the frequency of detection of the atom numbers, set for 
instance by varying the intensity of the detection light; 
from top to bottom we have NT/J = 1, 100, and 10000. 
The graphs on the left present the population imbalance 
as a function of time computed from Eq. ([7]), while the 
graphs on the right show the population imbalance z(t) 
versus the angle fit), Eq. ([5]), as points for all of the 5000 
time samples that were also used to plot z(t). 

First look at the middle data set with NT/J = 100. 
In the plot of z{t) all three behaviors anticipated for the 
classical system are apparent at once: symmetric oscilla- 
tions between the wells, and oscillations in which a ma- 
jority of the atoms stays in one or the other of the wells. 
A classical system starting from a random initial condi- 
tion and thereafter evolving deterministically would stick 
with one type of behavior. We attribute the switching to 
quantum noise and/or noise due to the back-action of 
the measurements. On the other hand, the plot of z{(p) 
is as if the system stayed close to the homoclinic orbits, 
something that the classical system is expected to do. 

The top data set differs in the frequency of observa- 



tions, NT/ J = 1. Wc qualitatively assign a period of the 
oscillations in population trapping, approximately half 
of the period of those oscillations in which the atoms go 
back and forth between the traps, as the characteristic 
time scale for the oscillations. Then the top panel cor- 
responds to about five photon counts per characteristic 
time, the middle panel to about five hundred. There is 
still a qualitative resemblance to the classical dynamics 
in the results for NT/ J ~ 1, but the quantitative agree- 
ment is quite poor compared to the NT/ J = 100 case. 
When the rate of observations is further decreased, the 
resemblance to classical dynamics deteriorates further. 
In a very real sense the observations cause the classical 
behavior: for the classical dynamics to emerge, the state 
has to be observed frequently enough that the classical 
behavior can be resolved. 

In the bottom data set we have NT/ J = 10"', which 
means that several photons are recorded per each atom 
during a characteristic time of the oscillations of the 
atoms between the traps. Although the noise is en- 
hanced, the classical behavior is still clearly discernible. 
Our interpretation is that, while increasingly aggressive 
measurements will eventually destroy the classical behav- 
ior, the system is highly resilient against the noise from 
the back-action of the measurements. 

Our particular model for continuous measurements of 
the atom number was picked for numerical expediency, as 
quantum trajectory simulations boil down to solving the 
Schrodinger equation. Real experiments on the instabil- 
ity would not necessarily comply. For instance, suppose 
one has inside a Mach-Zehnder interferometer light that 
is far-off resonance from any atomic transition, and that 
the interferometer is initially balanced so that all light 
comes out from one port. When one inserts an atomic 
sample inside one arm of the interferometer, the refrac- 
tive index from the atoms make light to come out of the 
formerly dark port. But the amplitude of this light is 
proportional to the number of the atoms, so that the in- 
tensity, and photon counting rate, is proportional to the 
square of the number of the atoms not atom number |l6| . 

The notion that measurements or other environmental 
influences may cause a quantum system behave approx- 
imately classically is by no means new [19[, and quan- 
tum trajectory type methods have been used to study 



the emergence of classical chaos in quantum systems |20l - 
|23| | . The unusual feature about our dynamically unstable 
model is that a dramatic macroscopic behavior is trig- 
gered entirely by measurement back-action. We do not 
enter into details of different measurement schemes, for 
one thing because we believe that within reason (many 
atoms, measurements not too weak or strong) classical 
physics always emerges |2l| - |23| . However, although the 
smaller atom numbers dictated by the constraints of the 
numerics make the results less striking, we have verified 
the same qualitative behaviors as in our explicit exam- 
ples also in the case when the photon counting rates are 
proportional to the squares of the atom numbers. 

Making the sides of the trap and the measurement 
scheme symmetric enough that quantum phenomena 
rather than technical imperfections dictate the behavior 
of the atoms may be a challenge in real experiments. On 
the other hand, the assumption that no scattered photon 
is lost without a detection event is clearly unimportant 
for the basic principle. For instance, the detection of the 
classical behavior suffers because of the less-than-unit ef- 
ficiency of a photon counter, but information about the 
atoms numbers gets conveyed to the detectors and there 
is a back-action from the measurements. Analysis of the 
effects of nonidealities of the system and of the measure- 
ments is an interesting problem area, but we leave it for 
future work. 

We have studied the behavior of a Bose-Einstein con- 
densate in a double- well trap under the conditions when 
the classical model of the system has an unstable steady 
state and exhibits large-scale nonlinear dynamics as a re- 
sult. A priori, the corresponding quantum system does 
no such thing, but when the state is monitored the clas- 
sical nonlinear dynamics nonetheless emerges. That a 
quantum mechanical calculation gives you nothing with- 
out a description (at least an implied description) of the 
measurements may be trivial in itself. We think, though, 
that we also have a practical point: When one deals 
with a system that is intermediate between quantum and 
classical mechanics, a careful description of the measure- 
ments is a necessary part of the analysis. 
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